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Background oriented schlieren (BOS) is a recent development of the schlieren and shadowgraph methods. 
The BOS technique has the ability to provide visualizations of the density gradient in both the axial and 
radial directions. The resultant magnitude of the density gradients allows for comparison with shadowgraph 
images. This paper first compares data obtained by the BOS and shadowgraph techniques at identical 
conditions in a free jet. The patterns and spacing of the shock trains obtained by the two techniques are 
found to be consistent with one another. This provides confidence in the shock spacing measurement by the 
BOS technique. Due to its simpler setup, BOS is then applied to investigate the shock spacing associated with 
the screech phenomenon, especially during stage jumps. Screech frequencies from a 37.6 mm convergent 
nozzle, as a function of jet Mach number (My), are shown to exhibit various stages. As many as eight stages 
are identified with the present nozzle over the range 1.0 < Mj <1.7. BOS images are acquired at various 
screech conditions and the shock spacing is examined as a function of Mj. 
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Nomenclature 

angle through which light is refracted 
refractive index gradient 
apparent pixel shift due to refraction 
jet Mach number 
design Mach number 


I. Introduction 

The aim of the present investigation is first to assess the effectiveness of measuring shock spacing using the BOS 
technique by comparing its results with data from the conventional shadowgraph technique and then applying it to 
study shock spacing associated with the screech phenomenon. The BOS technique and aspects of the screech 
phenomenon are reviewed in the following sections. 

A. Background Oriented Schlieren 

Background oriented schlieren (BOS) is a recent development of the schlieren and shadowgraph techniques, 
which are used to non-intrusively visualize density gradients. Schlieren and shadowgraph methods provide 
qualitative information of the first and second derivatives of the density in the flow, respectively. However, 
schlieren and shadowgraph techniques can be difficult, time consuming and costly, as they require large mirrors or 
lenses and precise alignment. Conversely, BOS captures the density field but only requires a CCD camera, light 
source, and a background, making it a much simpler and less costly option for flow visualization 1 . 

BOS is based on background distortion or an apparent movement of the background when imaged through a 
density field onto a detector plane 2 . Figure 1 illustrates the principle of BOS. An incoherent light source uniformly 
illuminates a background composed of a high-contrast random pattern that is imaged by a camera and lens system. 
Without any refractive index gradient (denoted ‘S’ in Fig. 1) in the optical path, light will travel from the 
background undisturbed. The dashed green line, along the optical axis in Fig. 1, illustrates an individual light ray 
traveling from a point on the background to the camera undisturbed in the case without a refractive index gradient 
present. When acquiring an image under this condition, the background will look undistorted. A refractive index 
gradient in the path causes light to refract. The solid green line in Fig. 1 shows the light ray being refracted by the 
gradient at ‘S’. Now when the background is imaged through the non-uniform density field, points on the 
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background will appear shifted by ‘Ax’ at the detector plane. Therefore, when acquiring BOS data it is necessary to 
take two images: a reference image prior to the refractive index gradient being present and an image after the light 
has refracted through the density field. For maximum sensitivity the refractive index gradient should occur halfway 
between the camera and the background 2 . 

The shift between the two acquired images caused by the refracted light can be calculated by correlation 
methods. Existing algorithms that have already been developed for particle image velocimetry (PIV) are used to 
determine the background image displacements within small sub-regions across the image 3 . Unique to BOS, the 
density gradient in two orthogonal directions in the plane of the background can be distinguished from one another. 


Light source 



Background 

Fig. 1 Illustration of the principle of BOS, showing light passing through a refractive index gradient (S) 
resulting in a displacement of Ax. This displacement appears as a pixel shift in the acquired images. 


B. Screech 

Supersonic jets operating at imperfectly expanded conditions generate noise, which is comprised of three basic 
components: turbulent mixing noise, broadband shock-associated noise and screech tones . 4 The latter two noise 
components are associated with the quasi -periodic shock cell structures formed in the jet plume. The screech 
phenomenon, involving discrete frequency tones, is the primary focus of the BOS application in this study. 

Alan Powell first reported detailed observations of screech tones and offered an explanation for their occurrence 
based on acoustic feedback phenomenon 5 . The feedback mechanism responsible for screech tone generation has 
continued to be studied for several decades 6,7 . The feedback loop begins in the thin jet mixing layer, near the nozzle 
lip, where it is most receptive to external excitation. Acoustic disturbances traveling upstream impinge on this area 
and excite the inherent instability waves. While propagating downstream the instability wave grows and when 
sufficiently amplified it interacts with one of the quasi-periodic shock cells, generating a disturbance. The 
disturbance travels outside of the jet back upstream, impinges on the jet mixing layer near the nozzle lip and excites 
new instability waves, thus closing the feedback loop. The process is accompanied by the generation of an intense 
tone that is referred to as screech. 

As the screech tone is emitted, the jet flow undergoes strong oscillations. The oscillation can be in toroidal mode 
or in helical/flapping mode. The flapping mode oscillations consist of equal amounts of left and right helical mode 
oscillations. It has been known that the dominant screech tone goes through mode-switching or stage-jumps as the 
Jet Mach number (Mj) is increased. The frequency within a given stage decreases continuously as Mj is increased. In 
previous studies with circular nozzles , 8,9 flow field oscillations in the axisymmetric (stages Ai and A 2 ), flapping 
(stages B and D), and helical (stage C) modes have been observed. However, the B and D stages may be unstable, 
switching between the flapping and helical oscillations. The staging behavior will become clearer with the frequency 
data presented in the following. However, it is fair to say that the screech mechanism has remained from being 
completely understood due to its unsteady, intermittent nature and its dependence on a variety of parameters. One of 
the unknowns is the behavior of the shock spacing during a stage jump. It is not clear if the shock spacing adjusts to 
accommodate the new wavelength associated with a new stage or if it remains unchanged while other parameters 
adjust to accommodate the new stage. In this paper, we examine the shock spacing variations using BOS. 
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II. 


Experimental Setup and Procedure 


The experiments were conducted in a large open jet facility at NASA Glenn Research Center, shown in Fig. 2(a). 
Compressed air is supplied through one end of the 760 mm diameter plenum chamber fitted with flow conditioning 
units. The air discharges through the nozzle on the other end of the plenum into the quiescent ambient of the test 
chamber. All data in this paper pertain to unheated flow. Figure 2(b) shows a 50.8 mm diameter convergent- 
divergent nozzle with a design Mach number (M D ) of 1.8 which was previously used to acquire shadowgraph data in 
connection with internal shock study. 10 BOS data obtained with this nozzle is first compared to the shadowgraph 
data. Figure 2(c) shows a 37.6 mm diameter convergent nozzle used for the subsequent screech study. 



Fig. 2 Experimental facility, (a) Open jet rig, (b) 50.8 mm, M D = 1.8, convergent-divergent nozzle used to 
compare BOS and shadowgraph data, (c) 37.6 mm convergent nozzle used to study screech phenomenon. 

A top-view diagram of the BOS setup with relevant dimensions is shown in Fig. 3. A 10 W, 530 nm light 
emitting diode (LED) is used as the incident light source for the system. Light from the LED is uniformly spread 
across the 0.3 m x 0.3 m background, consisting of a random speckle pattern of retro -reflective paint particles, 
which directs the incident light back onto itself and towards the camera. The retro -reflective paint consists of 
thousands of tiny glass beads bonded to the painted surface, creating a high-contrast background. The image of the 
illuminated background is focused onto the camera detector by a Nikon 200 mm lens. The scientific -grade CCD 
camera is a Princeton Instruments model EC 1 1000 which has a detector with 4008 x 2672 pixels and a pixel size of 
9 pm x 9 pm. The 12-bit images were captured with a 30 ms exposure time. The large detector size enabled capture 
of the desired field of view within the distances dictated by the constraints of the test facility. The background was 
placed at a distance of 1.74 m from the focusing lens of the camera. The distance from the camera lens to the nozzle 
was 0.984 m and the distance from the background to the nozzle was 0.756 m. 

To obtain the BOS data, image pairs were acquired over a range of Mj. A reference image was acquired prior to 
the flow being turned on and the data image was acquired upon reaching the desired condition. The pixel shifts, in 
both the axial (x) and radial (y) directions, between the two images were calculated using PIV software. Using 
Tecplot, the pixel shifts were plotted as contours, which enabled visualization of the density gradients in the axial 
and radial directions. The magnitude of the resultant density gradients was plotted and compared with the 
shadowgraph data taken at the same conditions. 
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Fig. 3 Top view schematic of the BOS setup. 


III. Results 


A. BOS Data Compared to Previously Acquired Shadowgraph Data 

As stated before, the BOS data images were acquired with the 50.8 mm diameter, M D = 1.8 nozzle (designated as 
‘M18L’ in Ref. 10). The images were obtained over the range 1 .2<M / <1 .7 at matching test conditions of the 
previous shadowgraph experiment reported in the cited reference. An example of the BOS data at Mj= 1.6 is shown 
in Fig. 4. The pixel shifts measured in the v-direction, corresponding to the axial density gradient, is shown in Fig. 
4(a). Pixel shifts in the y-direction, corresponding to the radial density gradient, are shown in Fig. 4(b). These 
figures illustrate the usefulness of BOS for obtaining individual components of the density field in a single snapshot. 
The magnitude of the two components is shown in Fig. 4(c); this is the quantity that was used for comparison with 
the shadowgraph data. 

Examples of shadowgraph images at M/=1.2, 1.4, and 1.6 are displayed in Figs. 5(a)-5(c), respectively. 
Corresponding BOS images, taken for the same values of Mj are displayed in Figs. 6(a)-6(c), respectively. It can be 
seen that the flow and shock boundaries are sharper in the shadowgraph data than in the BOS data. Also the BOS 
images contain some asymmetry. It is believed that with further refinement of the BOS technique, the image quality 
and asymmetry can be improved upon. Presently however, the overall agreement of the images obtained by the two 
techniques is good, particularly the inferred shock spacings as described further in the following paragraphs. 
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Fig. 4 Contour diagrams of BOS pixel shifts for the M D =1.8 nozzle at Mj- 1.6. (a) Shift in the axial (x) 
direction, (b) in radial (y) direction and (c) the resultant magnitude. 
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The shock spacing is examined in Fig. 7 for Mj = 1.4. The tips of the shock at the jet’s shear layer are marked by 
arrows. Average of the first two tip-to-tip spacings is taken to be the shock spacing. This is deemed justified because 
the spacing is not expected to change much within the first two cells, while an average provides a better 
representative value. As it can be seen from the labels in the figure, the average spacing obtained by the two 
techniques is in good agreement. 

It is important to note the subjectiveness in the determination of the shock tip location with either the BOS or the 
shadowgraph technique. The uncertainty in the shock spacing measurement with the BOS technique in Fig. 7 is 
estimated to be 7.6% whereas the uncertainty in the shadowgraph shock spacing is estimated to be 3.5%. The 
shadowgraph results in images that contain data at every pixel resulting in a spatial resolution of 0.248 mm/pixel, 
resulting in lower measurement uncertainty. Background oriented schlieren produces images that have an averaged 
data value for every 8 pixels by 8 pixels sub-region, resulting in a spatial resolution of 0.365 mm/subregion. Note 
that the camera used in the BOS experiment had much larger sensor with 10.5 times more pixels than the camera 
used in the shadowgraph experiment. Although, the spatial resolutions between the two techniques are comparable, 
the difference in the images can be attributed to the different setups, light sources, field of views, dynamic ranges 
and sensitivities to the density gradients present in these shock-laden flows. 

Figure 8 shows the comparison of shock spacing versus M ; over the range 1.2 <Mj <1.7 obtained by the two 
techniques for the M D = 1.8 nozzle. Shock spacing calculated from the shadowgraph and BOS images are consistent 
within the data uncertainty and follow the expected trend previously observed in Ref. 6. As Mj increases the shock 
spacing increases as well. These results provide confidence in the shock spacing measurement by the BOS 
technique. The results on the study of shock spacing in a screeching jet are discussed in the next section. 


Mi =1.4 



b.) 



Average shock spacing = 0.9 65D 
Fig. 7 Average shock spacing determined from (a) shadowgraph and (b) BOS for M D = 1.8 nozzle at M 7 =1.4. 



Fig. 8 Non-dimensional shock spacing calculated from shadowgraph 
and BOS images as a function of Mj for the M D = 1.8 nozzle. 
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B. Screech Data 

The screech phenomenon was studied with a 37.6 mm diameter circular convergent nozzle. The convergent 
nozzle was chosen since it pertains to most of the previous studies in the literature on the subject. In order to 
investigate the role of shock spacing in mode-switching or stage jumps, data were obtained with special attention to 
regions of stage jumps involving hysteresis. First, screech frequency data were acquired by spectrum analysis of the 
signal from a microphone positioned at 25° relative to the jet axis. Figure 9 shows the recorded screech frequencies 
plotted versus Mj depicting the various screech stages that were observed. Stages Ai and A 2 (axisymmetric), B and 
D (flapping), and C (helical) have been previously noted in the literature. Stages B and D have also been observed 
to be unstable and switch between the flapping and helical modes. The B stage (flapping) overlaps with C stage 
(helical) and the D stage (flapping) overlaps with C stage. In each of these overlap regions there is hysteresis, i.e., 
the stage jump occurs at a different location depending on whether Mj is increased or decreased. Thus, within a 
jump region, data could be obtained at exactly the same Mj but in different stages. Note that Fig. 9 also shows 
additional screech frequency stages not previously observed, to the best of the authors’ knowledge. These are 
labeled as B’, E, and F. Within the Arrange covered, previous studies noted only up to stage D. It appears the 
screech stages are dependent on details of the nozzle and the reason for this discrepancy remains unclear at this time. 

As stated before, when two screech stages overlap there is usually hysteresis. As Mj is increased the screech tone 
frequency in a particular stage decreases and jumps to a new stage at a particular Mj\ however, as Mj is subsequently 
decreased the return to the earlier stage takes place at a lower Mj. Examples of sound pressure level spectra are 
shown in Fig. 10; (the same color codes are used for the spectral traces as for the stages in Fig. 9). Figure 10(a) 
shows spectra for stages Ai and A 2 . For this stage jump hysteresis was not clear. On the other hand, Figs. 10(b)— (f) 
show spectra pairs in regions involving clear hysteresis for jumps in stages (B’, B), (B, C), (C, D), (D, E), and (E, 
F), respectively. Of note is Fig. 10(c) where stages B and C are captured at exactly the same frequency depending on 
increasing or decreasing Mj. 

It was theorized that the shock spacing would display an abrupt change at different stages of the Mj values where 
the hysteresis occurred. In order to investigate this, BOS data were acquired for conditions corresponding to the 
different data points of Fig. 9 and representative images are displayed in Figs. 11(a) - 11(1). Figure 11(a) shows an 
image without any perceptible screech tones present. Figures 11(b) - 11(f) show images for screech stages, Ai, A 2 , 
B’, B’+B, and B, respectively. The rest of the images in Fig. 11 show pairs in the overlapping (hysteresis) regions 
where the adjacent stages could be captured at the same Mj. Thus, Figs. 11(g)- 11(1) represent image pairs in the 
hysteresis regions for stages (C, D), (D, E), and (E, F), respectively. 



Mj 

Fig. 9 Fundamental screech frequency versus Mj for the 37.6 mm circular convergent nozzle. 
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M =1.085 M = 1.206 








Fig. 10 Sound pressure level spectra for various screech stages, (a) Spectra from the non-overlapping stages 
Ai and A 2 , (b) stages B’ and B, (c) B and C, (d) C and D, (e) D and E and (f) E and F. Suffix 6 U’ (or no suffix) 
in the legend denotes the spectra were recorded for increasing My, whereas the letter 6 D’ denotes the spectra 
were recorded for decreasing Mj. 
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Fig. 11 BOS images at different M j} (a) without any perceptible screech, (b) stage A l9 (c) stage A 2 , (d) stage B’, 
(e) stages B’+B, (f) stage B. (Figure continued on next page) 
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Figs. 11 BOS images at different My, (g,h) C and D stages within hysteresis loop for increasing M 7 and then 
decreasing Mj , respectively, (i,j), D and E stages within hysteresis loop, (k,l) E and F stages within hysteresis 
loop. 
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From the BOS images, the average shock spacing was determined in the manner previously described in section 
III. A. The resulting non-dimensional shock spacing measured as a function of Mj is shown in Fig. 12. It is apparent 
that the shock spacing follows a monotonic trend. As Mj increases, the shock spacing increases as well. 6 Large 
departures from this trend, beyond the uncertainty of the data, do not appear over the entire range of Mj covered in 
the experiment. At the beginning of the experiment, it was suspected that the shock spacing would display an abrupt 
change for different stages at the Mj values where the hysteresis occurred. No such change did occur. Therefore, it is 
inferred that the shock spacing is not the parameter that adjusts to accommodate a new frequency when a stage jump 
occurs. Other parameters, such as convection velocities in different segments of the feedback loop or the distance 
over which the instability wave develops nonlinearity, possibly adjust to provide for the stage jumps. 



Fig. 12 Non-dimensional shock spacing versus Mj 
corresponding to the frequency data shown in Fig. 9. 


IV. Conclusions 

Background oriented schlieren (BOS) images were compared with shadowgraph data acquired under the same 
conditions. The average shock spacing s obtained by the two techniques, in a convergent-divergent nozzle with 
design Mach number of 1.8, were shown to be consistent with one another. Although the resulting shadowgraph and 
BOS images had comparable spatial resolutions, there were differences in the images. The BOS images were not as 
crisp as the shadowgraph data and showed some asymmetry. This is thought to be due to lack of complete 
refinement of the BOS application in the present setup, e.g., background and the image processing. However, the 
shock spacing measured by BOS was found to be consistent with the shadowgraph data. This provided confidence in 
the shock spacing measurement by the BOS technique which was then used to investigate mode- switching behavior 
associated with the screech phenomenon. Screech frequencies from a circular convergent nozzle were collected 
over a jet Mach number range of 1.0 < Mj <1.7. The data exhibited as many as eight ‘stages’ within the Mj range 
covered. BOS images were acquired at different screech conditions with special attention to regions of stage jumps. 
The shock spacings from the BOS images were shown to follow a monotonic trend with varying Mj without any 
discontinuities across the stage jumps. Thus it was concluded that the stage jumps did not occur due to an 
adjustment of shock spacing. Other parameters, such as convection velocities in different segments of the feedback 
loop or the distance over which the instability wave develops nonlinearity, possibly adjust to account for the stage 
jumps. 
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